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ABSTRACT: Electroluminescence produced during avalanche development in gaseous avalanche 
detectors is an useful information for triggering, calorimetry and tracking in gaseous detectors. 
Noble gases present high electroluminescence yields, emitting mainly in the VUV region. The 
photons can provide signal readout if appropriate photosensors are used. 

Micropattern gaseous detectors are good candidates for signal amplification in high background 
and/or low rate experiments due to their high electroluminescence yields and radiopurity. In this 
work, the VUV light responses of the Gas Electron Multiplier and of the Micro-Hole & Strip Plate, 
working with pure xenon, are simulated and studied in detail using a new and versatile C++ toolkit. 
It is shown that the solid angle subtended by a photosensor placed below the microstructures de- 
pends on the operating conditions. The obtained absolute EL yields, determined for different gas 
pressures and as functions of the applied voltage, are compared with those determined experimen- 
tally. 



KEYWORDS: Detector modelling and simulations II (electric fields, charge transport, 
multiplication and induction, pulse formation, electron emission, etc); Ionization and excitation 
processes; Micropattern gaseous detectors (MSGC, GEM, THGEM, RETHGEM, MHSP, 
MICROPIC, MICROMEGAS, InGrid, etc); Scintillators, scintillation and light emission 
processes (solid, gas and liquid scintillators). 



* Corresponding author. 



Contents 



|I[ Introduction 

|2[ Simulation 

g GEM 

p.l| Simulation details 



3.2 VUV EL yield 



3.3 Solid angle 



3.4 Ratio between light and charge carriers [7] 



MHSP | 

p~T| Simulation details | 

VUV EL yield g 

[4.3| Spatial distribution of excitations [T(] 

Conclusions |l0| 



1. Introduction 

Cryogenic dual-phase liquid/gas as well as gaseous detectors have been or are currently being 
developed for direct detection of WIMPs (Dark Matter candidates) [Q, |3[ |j] and for neutrino 
less double beta decay (Ovj8j3) searches [||, [7j H|. These detectors consist of Time Projection 
Chambers (TPCs) operating with pure noble gases, where both primary VUV scintillation and 
primary ionization are produced by radiation interaction. The primary VUV light is used for start- 
of-event triggering while the electrons resulting from ionization are driven to a region where their 
signal is amplified. It is important to have the highest possible signal gain, given the low rate and/or 
high background nature of these experiments. Because of this, electroluminescence (EL) is usualy 
used for signal amplification, rather than charge multiplication. In this process, the electrons are 
accelerated by a suitable uniform electric field that, between collisions, supplies them with enough 
energy to excite but not to ionize atoms of the gas. For atmospheric pressure or higher, the excited 
atoms decay through the formation of excimers [g]. In the case of xenon, the result is the emission 
of VUV photons with a wavelength centered at 173 nm and with a Full-Width-at-Half-Maximum 
(FWHM) of 14 nm JTo|]. Further details about the electroluminescence process in pure noble gases 
can be found in reference [|T|] and references therein. 

The future ton-scale Dark Matter and Ovj3j3 experiments will introduce new challenges in the 



design of large- volume TPCs [12]. Namely, the increase in the number of the commonly used 
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PMTs for VUV light readout would lead to undesired high background levels. The high radiop- 
urity, simplicity, low cost and the possibility of constructing large areas, make of Micropattern 
Gaseous Detectors (MPGDs) good candidates for the amplification of the primary ionization sig- 
nal [ 13, 14 , 15 ]. In these devices, during avalanche multiplication, also VUV light is produced since 
both excitations and ionizations are induced by direct impact of electrons. It has been demonstrated 
that higher light gains than those obtained with parallel gaps can be achieved by using much lower 
voltages [16]. The VUV EL produced during the avalanches can be detected, e.g, by Geiger mode 
avalanche photodiodes (GM-APDs) |[l7|, |T8[ ], also characterized by a negligible natural radioactiv- 
ity ©. 

The simulation of EL produced in MPGDs is thus of major importance for the correct design 
of future detectors. In reference [ JIB] ] it was determined experimentally the absolute VUV EL yields 
of the Gas Electron Multiplier (GEM) and of the Micro-Hole & Strip Plate (MHSP). Recently, the 
yield of the thicker version of the GEM (the THGEM) was also measured [2C]. In this work, we 
use a recent developed and flexible Monte Carlo simulation toolkit [21] to estimate the properties 
of VUV EL produced in GEM and MHSP. The obtained EL yields are compared, for the case of 
pure xenon at different pressures, with the measurements. 



2. Simulation 



The toolkit used in this work is implemented in C++ and based on the Garfield [22] and Mag- 
boltz 8.9.7 [23, 24] programs. Through this flexible toolkit it is possible to follow the electrons 
through the gas in nearly arbitrary field geometries. Electrons are tracked at the microscopic level 
(atomic/molecular) using Monte Carlo procedures and cross -sections available in Magboltz. In the 
case of xenon, the gas studied in this work, beyond cross sections for elastic collisions and ion- 
izations, the program parametrizes the excited energy levels as functions of 50 energy groups. It 
is possible to obtain information about each excitation produced in the gas, namely: spatial po- 
sition, time, atomic level and spent energy. Additional details about the toolkit can be found at 
references [ pi] , 21]. 

We assume that every excited atom gives rise to the isotropic emission of a VUV photon [^, 
11, The toolkit and the model were previously validated for uniform electric field geometry in 
reference [21] by comparing the simulated results with measurements [^, ^7|] and earlier Monte 
Carlo work [25]. 

Electric fields created by microstructures like GEM and MHSP can not be determined analyt- 
ically due to their complex shape and to the presence of insulators. Thus, field maps were con- 
structed using ANSYS12® 



], a Finite Element Method (FEM) [g9j] program. The 3-D 10-Node 
Tetrahedral Electrostatic Solid finite element was used since it is suitable for modeling irregular 
meshes that need to adapt to curved edges. Since the hole and electrode patterns of the GEM and 
MHSP are periodic, we constructed unitary cells, being the longitudinal axis of the holes parallel 
to the zz' axis. Mirror symmetries applied in the xx' and yy' directions were used, creating maps of 
ideal microstructures with infinite areas. 



-2- 



3. GEM 



3.1 Simulation details 

For the simulations, a standard GEM was considered, consisting of a 50 /im thickness Kapton® foil 
covered with 5 jlim thickness copper layers on both sides and perforated with bi-conical holes of 50 
and 70 jUm diameter at the inner and outer apertures respectively, arranged in a regular hexagonal 
layout with an edge of 140 /im length. 

In reference [16], the primary charge was generated using a collimated (1 mm radius) 22.1 
keV X-ray beam and the VUV EL photons were measured by a Large Area Avalanche Photodi- 
ode (LAAPD) with an active diameter of 16 mm placed below the GEM. The drift and induction 
fields were 0.5 and -0.1 kVcm -1 respectively. The unusual reversed induction field was used ex- 
perimentally so that no charge readout below the GEM was needed and that the LAAPD could 
replace it. This experimental field setup was implemented in the simulations, as well as the beam 
characteristics. 

The primary electrons were released isotropically 250 /im above the GEM (measured in the 
zz' direction), in a region where the equipotential surfaces are completely flat and the electric field 
is uniform. The (x,y) coordinates of the starting point were sampled randomly according to the 
geometry of the X-ray beam, the corresponding absorption probability along the 8 mm length 
drift region and the transversal diffusion that the primary electrons undergo before arriving at the 
z = 250 plane. Each electron started with a kinetic energy sampled randomly according to 
the energy distribution given by Magboltz for the experimental drift field. A set of 3 ,000 primary 
electrons was simulated. Pure xenon at a temperature of 300 K and pressures of 1.0, 1.5 and 2.5 
bar was considered. 



3.2 VUV EL yield 

The calculated absolute VUV EL yield, i.e., the total number of VUV photons produced per pri- 
mary electron, is shown in Figure (lines) as a function of the potential applied between the top 
and the bottom electrodes of the GEM, Vgem- Results for the different simulated pressures are 
compared with those of the measurements reported in reference [|^] (open symbols). Specially 
for high voltages, the simulated VUV EL yield approaches the experimental results (which have 
an associated uncertainty of the order of 20 %). However, some differences between the data are 
present. 



3.3 Solid angle 

Not all of the VUV photons produced in the gas impinge the surface of the photodetector. The 
LAAPD has a subtended solid angle that need to be considered in the calibrations performed during 
the measurements. In reference [16] this solid angle was calculated considering that all the VUV 
photons were emitted from a plane corresponding to the bottom surface of the bottom electrode, 
being independent of the applied voltage and of the gas pressure. 

The plots of Figure || show the distribution of the z coordinate of the points where excitations 
occurred, for some applied potentials and for xenon at 1.0 bar and 300 K. In the same plots, it 
is included the same distribution for the ionization points. For a fixed pressure, at low voltages, 
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Figure 1 : Average number of VUV photons produced per primary electron as a function of the 
potential applied between the top and the bottom electrodes of the GEM, shown for different pres- 
sures of xenon at 300 K. The simulated results (lines) are compared with the total VUV EL yield 
reported in [|^] (open symbols). The latter were corrected for the actual fraction of solid angle (see 
Section 3.3 for further details) and the result is shown as closed symbols. 

the number of produced secondary charges is low (refer also to Figure |^) and the microstructure 
works almost only in EL mode, being ionizations rare. The regions where more atoms are excited 
reflect the volumes where higher electric fields exist. The fraction of excitations and thus, of VUV 
photons emitted in the upper half part of the holes is roughly the same as that in the lower half 
part. As the applied voltage increases, the size of the avalanche increases and a higher fraction of 
secondary electrons is created in its last stage (close to the bottom electrode), reflecting the cascade 
nature of avalanches. These secondary charges also excite atoms and thus induce the emission of 
VUV photons. The fraction of photons emitted in the upper half part of the holes is now much 
lower than that of photons originated in the lower half part. Since the GEM itself blocks a higher 
fraction of photons originated in the upper regions of the holes than of those from the lower region, 
the subtended solid angle increases with the applied voltage. 

For a given applied voltage, lower pressures mean higher gains, due to the increase of the en- 
ergy achieved by electrons between collisions. The ratio between the number of photons originated 
from the lower region of the holes and the number of photons coming from the upper region is, thus, 
higher for lower pressures. The subtended solid angle increases as the gas pressure decreases. 

In conclusion, the fraction of solid angle subtended by the LAAPD (or any other photosensor 
placed below the GEM), £2, is dependent on the gas pressure and on the applied voltage. Experi- 
mentaly it is difficult to determine Q. as function of these two quantities and thus, in reference Jl6| ] 
it was considered as being the same for all the conditions, giving more importance to photons 
produced in the lower region of the holes, which is appropriated only for high VUV EL yields. 

Using the information about the positions where excitations occurred, available through the 
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Figure 2: Distribution of the z coordinate of the points where excitations and ionizations occurred, 
shown for some voltages applied along the GEM's holes. For each case, a set of 10 primary 
electrons was simulated for xenon at 1.0 bar and 300 K. The edges of the microstructure are pro- 
portionally represented at each plot. The electrons enter the GEM hole for positive values of z and 
drift toward negative values. 



simulation toolkit, it was possible to estimate the solid angle for each of the simulated conditions, 
considering that a VUV photon is emitted isotropically from each excitation point. According 
to the detector geometry, the photon can impinge the photosensor surface, be blocked by the mi- 
crostructure or hit the detector walls. In the latter case no reflection is considered. The value of Q. is 
determined as being the ratio between the number of photons that impinge the photosensor surface 
and the total number of photons produced in the GEM. For ilustration, in Figure [| are represented 
the points where excitations occured during one avalanche happening for a gas pressure of 1.0 bar 
and VQgjyj = 500 V. The trajectories of the VUV photons that hit the LAAPD surface are also 



-5- 




-0.074-0.072 -0.07 -0.068-0.066-0.064-0.062 -0.06 
X [cm] 




0.018 0.02 0.022 0.024 0.026 0.028 0.03 0.032 
y[cm] 



Figure 3: Points where excitations occured during one avalanche happening in Xe at 1.0 bar and 
300 K and for Vq EM = 500 V. The trajectories of the VUV photons that hit the LAAPD surface 
are also included. The edges of the microstructure are proportionally represented at each plot. 
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Figure 4: Fraction of solid angle subtended by the LAAPD (full lines) as function of Vgem (left) 
and the absolute VUV yield (right). Curves for the different simulated pressures are shown. The 
plots also show the estimated values of the solid angle subtended by an infinit plane located below 
the GEM (dashed lines). 



Figure £la| shows the calculated Hasa function of Vgem, for the different considered pressures. 
The estimated values of the solid angle subtended by an infinit plane located below the GEM are 
also shown. These values represent the maximum light detection coverage that can be achieved. In 



Figure 4b, the values of D. are shown as functions of the light gain, i.e., the absolute VUV yield. 
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For higher gains, £1 starts to saturate since almost all the excitations happen at the exit of the hole, 



reason why Q, also approaches the value calculated in reference []16[], £2 = 0.32. 

The experimental values of the absolute VUV yield, corresponding in Figure [T] to the open 
symbols, were corrected for the actual fraction of solid angle. The corrected values are shown 
in the same figure as closed symbols. It is clearly visible the better agreement between corrected 
experimental results and Monte Carlo simulations, specially for lower voltages. The overall agree- 
ment is remarkable considering the uncertainty in the measured values (~20 %). The small re- 
maining differences can be attributed to the effect of the charging-up of the insulator surface, not 
considered in the simulations; to intrinsic errors in the electric field calculations, characteristic of 
the FEM method; to imperfections in the GEM geometry introduced during its construction; and 
to the uncertainty in the temperature at which the measurements were performed. However, the 
contribution of these effects seems to be small, given the achieved good agreement. 

3.4 Ratio between light and charge carriers 

The spatial distributions presented in Figure f| show that, in a GEM, the number of excitations 
produced during avalanches is much higher than the number of ionizations. The ratio between 
these two numbers is shown, as a function of Vqem and for the different simulated xenon pressures, 



in Figure |5a[ As reference, in Figure |5b| it is shown the average number of secondary electrons 
produced per primary electron, for the same conditions. 

For the same applied voltage, as higher pressures are considered, the maximum energy ac- 
quired by electrons between collisions decreases, since their free path is lower. Since the energy 
threshold for ionizations is higher than for excitations, the latter are more likely for higher pres- 
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Figure 5: a) Ratio between the number of excitations and the number of ionizations produced 
during avalanches developed in a GEM, shown as a function of the voltage applied between the top 
and the bottom electrodes, Vgem. and for different pressures of xenon at 300 K. b) Average number 
of secondary electrons produced per primary electron (including the latter) for the same conditions. 



-7- 



sures. Considering the same pressure and starting with low voltages (V = 200 V), the ratio initially 
grows since the electric field is increasing and the maximum kinetic energy reached by the electrons 
is enough for excitations but barely for ionizations. For a certain voltage, at some of the collisions, 
electrons can achieve an energy slightly higher than the ionization threshold. At these energies, the 
ionization cross section is higher than that for the sum of the excitations and thus the ratio starts to 
drop. 



4. MHSP 



4.1 Simulation details 

We simulated also the VUV EL response of a standard MHSP, consisting of a 50 /im thickness 
Kapton® foil covered with 5 /im thickness copper layers on both sides and perforated with bi- 
conical holes of 50 and 70 jltm diameter at the inner and outer apertures respectively, arranged in 
an hexagonal pattern. The hole-pitch was 200 /im along the xx' direction (perpendicular to the 
strips) and 140 /im along the yy' direction (parallel to the strips). The width of the cathodes and 
anodes was 100 and 20 jJ,m respectively. The shape of the microstructure is depicted in Figure ||. 




x [ 10" D m] 



Figure 6: Equipotential lines produced in the MHSP microstructure when Vct = 450 V and Vac = 
250 V are applied. Isolines are shown for potentials between 100 V and 800 V in steps of 10 V, 
according to the colorbar shown at the right. The edges of the microstructure are proportionally 
represented at the same plot. The view plane, which contains the page, is y = 35 jUm and contains 
the center of the represented hole. The magenta lines show the average paths travelled by electrons 
(drift lines) starting 250 above the microstructure, being the x coordinate equally distributed 
along an horizontal line that is contained in the view plane. Drift lines beginning at x < and 
at x > 200 jUm are focused into the holes which centers are contained in a plane 70 /im behind the 
view plane. 
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The detector geometry of reference jl6| ] was considered. There, the primary charge was extracted 
from a semi-transparent Csl photocathode by photons emitted from a UV lamp. The VUV EL 
photons produced during avalanches extracted photoelectrons from a reflective Csl photocathode. 
The electrons were then collected by a metal mesh, being the induced current proportional to the 
number of produced VUV photons. The drift and induction fields were 0.1 and -0.1 kVcm -1 
respectively. 

The primary electrons were released isotropically 250 above the MHSP (measured in 
the zz! direction), in a region where the equipotential surfaces are completely flat and the electric 
field is uniform. The (x,y) coordinates of the starting point were sampled randomly according to 
an uniform distribution along the active area of the microstructure. Each electron started with a 
kinetic energy sampled randomly according to the energy distribution given by Magboltz for the 
drift field. A set of 500 primary electrons was simulated. Pure xenon at 1.0 bar and 300 K was 
considered. 

4.2 VUV EL yield 

The calculated total VUV EL yield is shown in Figure |as a function of the total voltage, Vtotai- 
This quantity corresponds to the sum of the potential applied between the cathode and the top 
electrodes, Vct, with the potential applied between the anode and the cathode, Vac- The results are 
compared with the measurements. The simulated and the experimental yields obtained for the case 
of GEM are also included in the plot. The two exponential regimes typical of the MHSP are clearly 
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Figure 7: Average number of VUV photons produced per primary electron in a MHSP as a function 
of the total applied voltage (Vtotai = Vct + Vac)- The results are shown for xenon at 1.0 bar and 
300 K. Those obtained through the simulation (green line) are compared with the total VUV EL 



yield reported in [16] (filled squares). The results of simulations and measurements (uncorrected 
for the variation of £2) for the case of GEM are also included in the plot. 
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visible in the Monte Carlo curve: the first corresponds to the increase on Vct while Vac = V, and 
the second to an increasing Vac while the value of Vct is kept constant. 

The Monte Carlo absolute VUV EL yield doesn't agree with the measured one. However, 
it is close to the experimental and simulated GEM's yields for the points where Vac = V (the 
initial exponential regime). Both the GEM and the MHSP considered in this work have exactly the 
same hole shape and dimensions, having the induction fields in both cases the same orientation. 
Thus, for the case of V A c = V, the VUV EL yields of the MHSP and the GEM should be very 
similar, which is not the case for the measurements. Experimentaly, for the case of the MHSP, 
the gas purification was achieved by diffusion, unlike as in the case of the GEM, for which it was 
achieved by convection. The xenon purity was thus affected, being possible that the measurement 
underestimated the VUV EL yield of the MHSP. It is known that small amounts of impurities have 
a big quenching effect in the production of VUV EL ligh in pure noble gases [|26|]. 

For the second regime (Vac > V), the experimental yield is now higher than the simulated 
one. This suggests imperfections in the construction of the microstructure, which is of peculiar 
difficulty. Namely, variations in the shape and dimensions of the anode have an important impact 
on the overall gain, both in terms of charge and light. Figure || shows the equipotential lines 
produced in the MHSP when Vct = 450 V and Vac = 250 V are applied, as well as the drift lines 
of some electrons. As shown, the highest electric fields exist near the anode, where the density of 
drift lines is also higher. Here, the drifting electrons have the highest probability of producing new 
ionizations. Besides that, this region consists of the second stage of charge amplification and thus, 
small variations in the field shape or amplitude have big effects on the overall gain. 

4.3 Spatial distribution of excitations 

Figure || shows the distribution of the (x, z) coordinates of the points where excitations occurred, 
for Vct = 450 V and Vac = 250 V (the maximum voltages considered in the simulations). The two 
regions with higher intensity reflect the two charge multiplication stages, one in the exit of the hole 
and the other in the vicinity of the anode strips. The presence of these two amplification stages is 
one of the figures of merit of the MHSP microstructure. 

5. Conclusions 

In this work, we estimated the total VUV EL yield of the GEM and MHSP microstructures working 
in pure xenon, through simulations at the microscopic level performed using a flexible C++ toolkit 
based on the Garfield and Magboltz programs. 

We demonstrated that the solid angle subtended by a photosensor placed below the GEM 
depends on the operating conditions. For the same pressure, as the voltage applied across the holes 
increases, the solid angle increases. For a fixed voltage, as the xenon pressure increases, the solid 
angle decreases. After correction for the solid angle variation, we were able to approximately 
reproduce the experimental measurements. 

The results obtained for the case of the MHSP are not consistent with measurements, although 
the typical two exponential amplification regimes are present in the Monte Carlo results. However, 
the comparison between the experimental yields for GEM and MHSP, when the voltage between 
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Figure 8: Distribution of the (x,z) coordinates of the points where excitations occurred, shown for 
V CT = 450 V and Vac = 250 V. The distribution refers to 2 primary electrons undergoing multipli- 
cation in a MHSP working with xenon at 1.0 bar and 300 K. The edges of the microstructure are 
proportionally represented at the same plot. 



the anode and the cathode of the latter is null, are unexpectedly very different. In these condi- 
tions, since the hole and induction field geometries are very similar, the MHSP should give similar 
response to that of the GEM. The disagreement suggests an experimental underestimation of the 
yield, supported by the used gas purification method: diffusion rather circulation. 
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